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vertical  columns  (or  piers)  and  horizontal  beams  (or  spandrels).  However,  with 
special  modeling  techniques,  very  complex  situations  may  be  considered.  A 
special  shear  panel  element  is  developed  to  enable  modeling  of  discontinuous 
shear  walls  and  shear  walls  with  arbitrary  openings.  A  diagonal  bracing  system 
to  model  braced  frames  (X-braced,  K-braced,  or  eccentrically  braced  systems)  is 
also  presented. 

The  column,  shear  panel,  and  diagonal  formulations  include  the  effects  of 
bending,  axial,  and  shear  deformations.  Bending  and  shear  deformations  are  also 
included  in  the  beam  formulation;  however,  the  effects  of  axial  deformations  are 
neglected. 

The  effects  of  the  finite  dimensions  of  the  beams  and  columns  on  the 
stiffness  of  a  frame  or  shear  wall  system  are  automatically  included. 

The  buildings  may  be  unsymmetrical  and  nonrectangular  in  plan.  Torsional 
behavior  and  interstory  compatibility  are  accurately  reflected  in  the  results. 

Four  independent  vertical  and  two  independent  lateral  static  load  condi-'* 
tions  are  possible  in  any  one  run.  These  six  static  load  conditions  may  be 
combined  in  any  ratio  to  each  other  or  to  a  lateral  dynamic  earthquake  input 
which  may  be  specified  as  a  time-dependent  ground  acceleration  or  as  an 
acceleration  response  spectrum. 

Three-dimensional  mode  shapes  and  frequencies  are  evaluated. 

The  unique  solution  procedure  used  by  CTABS80  considers  the  frame  and 
shear  walls  as  substructures,  reduced  with  a  modified  wave  front  technique. 

This  method  results  in  a  significant  reduction  in  the  program  data  preparation, 
computational  effort,  and  storage  requirements. 

The  consecutive  levels  of  each  of  the  individual  frames  can  be  arbi¬ 
trarily  connected  to  any  (sequential  by  not  necessarily  consecutive)  level 
of  the  structure,  thereby  making  it  possible  for  frames  tc  bypass  certain 
story  levels.  This  option  gives  the  program  the  capability  tc  model  partial 
diaphragms  and  multidiaphragms  at  any  level. 

A  user's  guide  for  the  program  is  presented  in  Waterways  Experiment 
Station  Instruction  Report  K-81-9. 


This  report  presents  the  theoretical  basis  for  a  computer  program 
called  CTABS80  that  can  be  used  for  static  and  dynamic  analysis  of  multi¬ 
story  frame  and  shear  wall  buildings.  Dr.  E.  H.  Wilson,  University  of 
California,  Berkeley,  was  responsible  for  developing  the  original  version  of 
the  program  (TABS),  sponsored  mainly  by  a  National  Science  Foundation  Grant. 

Modifications  to  the  program  to  make  it  a  more  useful  tool  for  Corps 
of  Engineers'  personnel  were  made  by  Mr.  Ashraf  Habibullah,  Computers/ 
Structures  International,  Oakland,  Calif.  His  work  was  sponsored  with 
funds  provided  to  the  Automatic  Data  Processing  (ADP)  Center,  U.  S.  Army 
Engineer  Waterways  Experiment  Station  (WES),  Vicksburg,  Miss.,  by  the 
Military  Programs  Directorate  of  the  Office,  Chief  of  Engineers,  U.  S.  Army 
(OCE),  under  the  Computer-Aided  Structural  Engineering  (CASE)  Project. 

This  report  and  a  companion  user's  guide  for  CTABS80  are  the  work  of 
Dr.  Wilson  and  Messrs.  H.  H.  Dovey  and  Habibullah. 

Specifications  for  the  modifications  to  TABS  were  provided  by  the 
members  of  the  CASE  Task  Group  on  Building  Systems.  The  following  were 
members  of  the  Task  Group  (though  all  may  not  have  served  for  the  entire 
period)  during  the  period  of  modifications  to  the  program: 

Mr.  Dan  Reynolds,  Sacramento  District  (Chairman) 

Mr.  Jerry  Foster,  Baltimore  District 
Mr.  Joseph  Hartman,  St.  Louis  District 
Mr.  David  Illias,  Portland  District 
Mr.  Sefton  Lucas,  Memphis  District 
Mr.  Jun  Ouchi,  Pacific  Ocean  Division 
Mr.  David  Raisanen,  North  Pacific  Division 
Mr.  Pete  Rossbach,  Baltimore  District 
Mr.  James  Simmons,  Baltimore  District 
Mr.  Ollie  Werner,  Middle  East  Division 
Mr.  Gene  Wyatt,  Mobile  District 

Dr.  N.  Radhakrishnan,  Special  Technical  Assistant,  ADP  Center,  WES, 
and  CASE  Project  Manager,  and  Mr.  Paul  K.  Senter,  Computer-Aided  Design 
Group  (CADG),  ADP  Center,  coordinated  and  monitored  the  work.  Ms.  Deborah 
K.  Martin,  CADG,  supported  the  Task  Group  in  changing  the  program  to  accept 
free-field  input.  Mr.  Seymour  Schneider,  Military  Programs  Directorate, 
was  the  OCE  point  of  contact.  Mr.  Donald  L.  Neumann  was  Chief,  ADP  Center. 

Directors  of  WES  during  this  period  were  COL  J.  L.  Cannon,  CE,  COL 
N.  P.  Conover,  CE,  and  COL  T.  C.  Creel,  CE.  Technical  Director  was  Mr. 

F.  R.  Brown. 
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CONVERSION  FACTORS,  INCH-POUND  TO  METRIC  (SI) 
UNITS  OF  MEASUREMENT 


i 


Inch-pound  units  of  measurement  used  in  this  report  can  be  converted  to 
metric  (SI)  units  as  follows: 


Mul tiply 

By 

To  Obtain 

feet 

0.3048 

metres 

inches 

2.54 

centimetres 

kips  (1000  lb  force) 

4.44 8222 

kilonewtons 

kips  (force)  per  foot 

14.593904 

kilonewtons  per  metre 

pounds  (force)  per  square  foot 

47.880263 

pascals 

pounds  (mass)  per  cubic  foot 

16.01846 

ki lograms  per  cubic 
metre 
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THEORETICAL  BASIS  FOR  CTABS80:  A  COMPUTER 
PROGRAM  FOR  THREE-DIMENSIONAL  ANALYSIS 
OF  BUILDING  SYSTEMS 


CHAPTER  I:  INTRODUCTION 


A.  Purpose 

This  report  presents  the  theoretical  basis  for  CTABS80,  a  computer  program 
for  the  linear  three-dimensional  structural  analysis  of  multistory  frame 
and  shear  wall  buildings  subjected  to  static  and  dynamic  loadings.  A  user's 
guide  for  the  program  is  presented  in  Waterways  Experiment  Station  (WES) 
Instruction  Report  K-81-9^8^. 

B.  General-Purpose  Programs  for  Structural  Analysis 

There  are  many  two-  and  three-dimensional  computer  programs  for  the  linear 

(1  2) 

analysis  of  complex  structures  .  Most  of  these  programs  can  be  used 

for  the  static  and  dynamic  analysis  of  multistory  frame  and  shear  wall 
buildings.  However,  most  of  these  programs  do  not  give  special  recognition 
to  the  fact  that  building  systems  in  themselves  are  a  very  special  class  of 
structures  from  an  analytical  point  of  view.  The  following  are  some  of  the 
characteristics  that  are  inherent  in  the  nature  of  a  building  analysis  that 
a  general-purpose  analysis  program  may  not  recognize,  thereby  resulting  in 
significant  losses  in  man-hours,  computer  time,  and  possibly  accuracy: 

1.  Most  buildings  are  of  simple  geometry  with  horizontal  beams  and 

vertical  columns.  A  simple  rectangular  grid  can  define  such  a  geometry 


4 


I 


tj 


vertical  columns.  A  simple  rectangular  grid  can  define  such  a  geometry 
with  minimal  input.  See  Figure  1. 

2.  Many  of  the  frames  and  shear  walls  are  typical.  Most  general-purpose 
programs  do  not  recognize  this  fact;  therefore,  the  input  may  be  large, 
and  some  internal  calculations  may  be  unnecessari ly  duplicated. 

3.  The  in-plane  stiffnesses  of  the  floor  systems  of  most  buildings  are 
very  high.  General-purpose  programs  do  not  necessarily  recognize  this, 
resulting  in  a  set  of  equilibrium  equations  which  may  be  very  large, 
and  thereby  causing  an  increase  in  computation  effort  by  a  factor  of  10 
to  100.  Also,  numerical  errors  may  be  introduced  since  the  in-plane 
floor  stiffnesses  are  several  orders  of  magnitude  greater  than  the 
story- to-story  stiffnesses  of  the  structure.  Since  these  two  stiffnes¬ 
ses  are  added  in  a  direct  stiffness  approach,  double  precision  may  be 
required  in  the  solution. 

4.  The  loading  in  building  systems  is  of  a  restricted  form.  Loads,  in 
general,  are  either  vertically  down  (dead  or  live)  or  lateral  (wind  or 
seismic).  The  vertical  loads  are  usually  applied  on  the  beams,  and  the 
lateral  loads  are  generated  at  the  floor  levels. 

5.  In  many  buildings,  the  dimensions  of  the  members  are  large  and  have  a 
significant  effect  on  the  stiffness  of  the  frame.  Therefore,  correc¬ 
tions  need  to  be  applied  to  the  member  stiffnesses.  Most  general- 
purpose  programs  work  on  center- line  dimensions,  and  stiffness 
corrections  are  usually  very  tedious  to  implement. 

6.  In  the  dynamic  analysis  of  buildings,  the  mass  of  the  structure  can  be 
accurately  lumped  at  the  floor  levels.  Recognizing  this  fact 
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significantly  reduces  the  size  of  the  eigenvalue  problem  to  be  solved. 


7.  Various  code  loading  requirements  necessitate  special  options  that 
allow  convenient  combinations  of  the  vertical  and  lateral  static  and 
dynamic  loadings.  Also,  the  member  forces  need  to  be  printed  out  at 
the  support  faces  of  the  members.  Such  transformations  are  not 
automatic  in  general-purpose  programs. 

8.  It  is  desirable  to  have  a  building  analysis  computer  output  printed  in 
a  special  format;  i.e.,  in  terms  of  a  particular  frame,  story,  column, 
and  beam.  Also,  special  output  such  as  story  shears  may  be  desirable. 

In  light  of  the  above-mentioned  and  other  reasons,  the  need  for  special- 
purpose  programs  for  building  analysis  is  apparent. 

C.  Special-Purpose  Programs  for  Building  Analysis 

Various  programs  have  been  developed  at  the  University  of  California  at 

Berkeley  for  the  linear  analysis  of  multistory  buildings  in  the  past  two 
(4  5  6l 

decades  '  ’  ’  .  These  programs  have  been  used  in  the  profession  on  many 

major  structures  in  many  different  countries.  One  of  the  major  reasons  for 

(I  2  3) 

the  development  of  computer  program  TABS  was  the  direct  "feedback" 

from  the  profession  in  the  use  of  these  programs. 

The  first  of  these  programs,  FRMSTC,  is  a  static  load  analysis  program  for 
symmetrical  buildings  with  parallel  frames  and  shear  walls.  Lateral  mode 
shapes  and  frequencies  are  also  evaluated. 

Program  FRMDYN  is  the  same  as  FRMSIC  except  that  the  load  input  is  ground 
accelerations  due  to  a  specified  earthquake.  Time-dependent  displacements 
and  member  forces  are  produced  but  are  not  combined  with  static  loads. 
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Program  LATERAL  is  an  extension  of  FRMSTC  to  the  static  analysis  of  a  sys¬ 
tem  of  frames  and  shear  walls  which  are  not  parallel.  Three  degrees  of 
freedom  exist  at  each  story  level.  This  program  does  not  have  dynamic 
options . 

The  first  version  of  TABS  was  released  in  1972,  with  the  intent  of  replac¬ 
ing  the  computer  programs  described  above.  CTABS80  is  an  enhanced  version 
of  the  original  version  of  TABS  and  is  intended  to  supercede  other  enhanced 
versions  such  as  XTABS  and  TABS77. 

The  computer  program  ETABS  was  released  in  1975.  The  program  allows 
three-dimensional  frame  input  in  which  common  column  compatibility  is  en¬ 
forced.  The  input  data  are  more  complex  than  those  of  TABS,  and  use  of 
this  program  is  only  recommended  if  common  column  compatibility  is 
important. 

For  buildings  with  other  complexities,  such  as  discontinuous  or  flexible 
diaphragms,  sloped  diaphragm,  nonrectangular  framing  systems,  etc.,  a 
general -purpose  program  such  as  SAPIV  or  EASE2  is  still  the  most 
appropriate  solution  tool. 

D.  Disclaimer 

Considerable  time,  effort,  and  expense  have  gone  into  the  development  and 
documentation  of  CTABSdO,  and  the  program  has  been  thoroughly  tested  and 
used.  In  using  the  program,  however,  the  user  accepts  and  understands  that 
no  warrantly  is  expressed  or  implied,  either  by  the  sponsors,  the  develop¬ 
ers,  or  the  distributors,  as  to  the  accuracy  or  the  reliability  of  the  pro¬ 
gram.  The  user  must  clearly  understand  the  basic  assumptions  of  the 
program  and  must  verify  his  own  results. 


8 


CHAPTER  II:  STRUCTURAL  IDEALIZATION 


An  exact  three-dimensional  structural  analysis  is  required  for  only  a 
limited  number  of  buildings.  For  the  majority  of  buildings  the  following 
approximations  can  be  made.  These  approximations  greatly  simplify  the 
preparation  of  input  data  and  significantly  reduce  the  computational 
efforts  associated  with  the  analysis  of  the  structure. 

1.  The  structure  is  idealized  as  an  assemblage  of  vertical  planar 
"frames1'.  A  frame  consists  of  m  columns  and  (m-1)  beams.  As  long 
as  shear  and  bending  deformations  are  included  in  all  members  there 
is  no  need  to  distinguish  between  a  shear  wal 1/spandrel  versus  a 
beam/column  system.  See  Figure  2. 

2.  The  out-of-plane  stiffnesses  of  all  frames  are  assumed  to  be  zero. 
Therefore  each  column  at  every  floor  has  two  degrees  of  freedom,  a 
vertical  displacement  and  a  rotation.  In  addition,  there  is  one 
lateral  degree  of  freedom  at  every  floor  level  of  the  frame. 

3.  Each  floor  is  modeled  as  a  horizontal  diaphragm.  This  diaphragm  is 
assumed  to  be  infinitely  stiff  in-plane.  The  out-of-plane  stiffness 
of  this  diaphragm  is  neglected.  Bending  stiffness  of  the  floors  may 
be  included  approximately  in  the  modeling  of  the  individual  frames. 

It  is  apparent  that  axial  deformation  is  not  permitted  in  the  beams. 
Floor  levels  must  be  the  same  for  all  frames.  See  Figure  2. 

4.  The  floor  diaphragm  connects  all  the  frames  together  at  the  corre¬ 
sponding  level.  The  connection  is  only  in  a  lateral  sense.  The 
frames  otherwise  are  completely  independent  of  each  other.  This 
also  means  that  compatibility  is  not  enforced  with  regard  to 
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Figure  2.  Typical  story  level 
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displacements  at  columns  which  are  common  to  more  than  one  frame. 

Thus  axial  deformations  in  common  columns  will  not  be  the  same.  As 
for  joint  rotations,  if  the  frames  with  common  members  are  perpendic¬ 
ular  in  plan  view,  then  the  rotations  are  uncoupled.  This  assumption 
invalidates  the  program  for  use  in  the  analysis  of  structures  in 
which  the  tubular  effect  or  common  column  compatibility  is  important. 

5.  Vertical  loads  are  applied  to  each  frame  on  a  tributary  area  basis. 
The  diaphragm  will  not  transfer  any  vertical  load  from  one  frame  to 
another.  However,  no  frame  can  sidesway  independently  without 
engaging  the  other  frames. 

6.  Lateral  loads  are  applied  as  loads  for  the  complete  floor  at  each 
level.  The  loads  are  applied  at  specified  locations  on  the  floor 
diaphragms  and  get  distributed  to  the  various  frames  in  accordance 
with  their  corresponding  stiffnesses  and  locations. 

A.  The  Frame  Substructure 

The  elevation  of  a  typical  frame  is  shown  in  Figure  3  .  The  frame  is 
basically  of  rectangular  geometry  with  vertical  column  center  lines  and 
horizontal  floor  levels  as  the  basic  reference  lines  for  the  description 
of  the  frame. 

The  frame  is  an  assemblage  of  column,  beam,  panel,  and  bracing  elements. 
Vertical  loading  is  applied  to  the  individual  frames  by  means  of  loading 
patterns  associated  with  each  beam. 

The  coljmn  and  beam  elements  have  options  for  rigid  offsets  at  each  end 
to  compensate  for  the  effects  of  the  finite  dimensions  of  the  members  on 


11 


Bea 


the  stiffness  of  the  system.  The  procedure  used  to  set  the  lengths  of 
these  rigid  offsets  is  presented  later  in  this  report. 


(i).  Individual  Member  Stiffnesses 

The  complete  stiffness  matrix  of  each  frame  is  assembled  by  the  direct 
stiffness  technique.  This  involves  calculating  the  local  stiffness  matrix, 
k,  for  each  member  along  with  a  transformation  matrix,  a,  which  transforms 
the  local  displacements  and  forces,  S,  to  global  displacements  and 
forces,  r,  R 

or:  <t>  =  a  r  also:  S  =  k  <j> 

S  =  a  R  and:  R  =  K  r 

where  K  is  the  stiffness  matrix  in  global  coordinates. 

Substituting  <f>  =  a  r  and  S  =  a  R  into  S  =  k  4>  we  get: 

a  R  =  k  a  r 

Premultiply  both  sides  by  a^  and  recognizing  a^  a  =  I  we  get: 

aT  a  R  =  aT  k  a  r 

R  =  a^  k  a  r 

As  R  =  k  r  we  get  K  =  aT  k  a 

Thus  knowing  the  local  stiffness  matrix  k  and  the  coordinate  transformation 

matrix  a  the  global  stiffness  matrix  may  be  evaluated. 

The  a  and  k  matrices  for  the  column,  beam,  panel,  and  brace  elements  are 
presented  in  Figures  4,  5,  6,  and  7,  respecti vely. 
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The  column  element  formulation  accounts  for  bending,  axial  and  shear  defor¬ 
mations.  The  basic  stiffness  matrix  for  such  an  element  is  shown  in 
Figure  8.  The  column  ends  have  options  for  rigid  offsets.  Figure  4  shows 
the  six  degrees  of  freedom  associated  with  the  column  element  and  the 
deformation  displacement  transformation  matrix  linking  the  frame  joint 
displacements  to  the  column  end  deformations. 

The  beam  element  formulation  is  similar  to  that  of  the  column  except  that 
the  axial  force  component  is  dropped  leaving  a  stiffness  and  transformation 
matrix  as  shown  in  Figure  5. 

The  panel  element  formulation  is  basically  the  same  as  that  of  the 
column  except  that  each  rotational  degree  of  freedom  is  transformed  into 
the  two  vertical  displacements  of  the  column  lines  bounding  the  panel 
element  at  each  corresponding  level.  The  axial  degrees  of  freedom  of 
the  panel  are  also  transformed  as  being  an  average  of  the  vertical 
degrees  of  freedom  of  the  two  column  lines  bounding  the  panel  element  at 
each  corresponding  level.  The  degrees  of  freedom  of  the  frame  associated 
with  the  panel  are  therefore  all  translational .  No  rigid  offsets  are  used, 
see  Figure  6. 

The  diagonal  element  formulation  is  exactly  the  same  as  that  of  the 
column  except  that  the  brace  is  inclined  and  no  rigid  offsets  are  used. 

See  Figure  7. 
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The  complete  stiffness  matrix  for  each  frame  has  two  degrees  of  freedom 
for  each  beam-column  intersection  and  one  lateral  degree  of  freedom  per 
story. 

(ii).  Lateral  Frame  Stiffness 

With  the  frame  degrees  of  freedom  appropriately  ordered,  the  frame  equi¬ 
librium  equations  have  the  form  shown  in  Figure  9  .  Where  N  is  the 
number  of  stories  in  the  frame,  r^  is  the  vector  of  joint  displacements 
(that  is  vertical  displacement  and  rotation)  at  story  level  n  and  r^  is 
the  vector  of  lateral  story  displacements.  Lateral  loads  are  applied  to 
the  complete  structure  and  are  considered  when  the  lateral  stiffness 
matrix  for  the  complete  building  is  assembled.  Gaussian  elimination  is 
performed  on  the  full  system  up  to  and  including  the  equations: 

-N  =  ?n_i  Dn-1  +  -N  -N  +  -N  -L 
The  last  N  equations  (r^  is  a  vector  of  order  N)  may  now  be  written  as: 

?L  =  -L  rL 

The  vector  is  the  lateral  load  submatrix  of  the  frame  and  is 
modified  by  the  elimination  process  due  to  vertical  loading  on  the  frame. 
These  terms  represent  the  sidesway  effects  under  vertical  loading.  The 
matrix  clearly  represents  the  frame  lateral  stiffness  matrix;  i.e., 
the  stiffness  matrix  of  the  frame  in  terms  of  only  the  lateral  story 
displacements. 

Within  the  computer  program  the  following  approach  is  adopted  in  order  to 
reduce  storage  requirements.  The  assembly  and  reduction  process  is 
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Figure  9.  Complete  equation  system  of  frame  substructure 
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.... 


carried  out  systematically  story  by  story  from  the  top  of  the  structure 
such  that  at  any  level,  n,  we  consider  the  system  shown  below: 
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where  the  prime  indicates  that  the  submatrices  may  have  been  modified 
by  previous  elimination. 


At  each  level  the  following  steps  are  performed: 

a.  Add  in  the  individual  member  stiffnesses  for  level  n. 
These  are  shown  below: 
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b.  Perform  the  elimination  on  the  equations  of  the  uppermost  partition 
in  the  equations  above 

c.  Save  these  reduced  equations  for  subsequent  back-substitution. 

d.  Rearrange  the  submatrices  in  the  equation  above  appropriately  in 
order  to  proceed  to  the  next  level.  This  rearrangement  is  as 
follows: 
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e.  Repeat  the  above  steps  for  the  next  level.  Thus  after  the  elimination 
is  completed  for  joint  displacements  at  all  story  levels,  we  are  left 
with  the  lateral  stiffness  matrix  for  the  frame. 

(iii).  Rigid  Joint  Offset  For  Beams  and  Columns 

The  deformations  within  the  joint,  an  area  bounded  by  the  finite 
dimensions  of  any  beam  and  column  intersection  (shown  shaded  in  Figure  10) 
are  neglected.  In  other  words,  this  area  is  assumed  to  be  an  infinitely 
rigid  rectangular  diaphragm. 

This  is  achieved  by  providing  rigid  offsets  at  the  ends  of  the  beams  equal 
in  length  to  one  half  of  the  widths  of  the  column  below  at  each  correspond¬ 
ing  end.  Rigid  offsets  are  also  provided  at  each  end  of  the  columns  equal 
to  the  depth  of  the  larger  of  the  beams  on  either  side  of  the  column  at 
the  corresponding  level. 

It  has  been  found  that,  in  general,  a  reduction  in  the  lengths  of  the 
rigid  offsets  to  compensate  for  some  deformation  that  may  exist  in  the 
joint  is  justifiable  and  gives  better  results,  especially  in  cases  where 
the  member  dimensions  are  substantial.  See  Reference  9. 

Reduction  of  the  rigid  link  dimension  has  been  coupled  to  the  size  of  the 
member.  In  other  words,  the  rigid  link  is  calculated  as  described  above 
and  then  is  reduced  by  25%  of  the  dimension  of  the  member,  at  each  end. 
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1  lustration  of  rigid  zone  reduction  methodology 


Thus  the  beam  rigid  links  are  reduced  by  25%  of  the  beam  depth  and  the 
column  rigid  links  are  reduced  by  25%  of  the  column  width  at  each  end. 
The  reduction  cannot,  of  course,  result  in  a  negative  rigid  link  length. 
This  reduction  procedure  is  optional.  If  no  column  widths  or  beam 
lengths  are  input  the  rigid  link  lengths  degenerate  to  zero  and  the 
analysis  is  carried  out  on  the  frame  grid  line  basis. 


B.  The  Complete  Structure 

In  order  to  combine  the  frame  lateral  stiffness  matrices  into  a  complete 
structure  lateral  stiffness  matrix,  each  of  the  frame  stiffnesses  must 
be  transformed  to  a  common  displacement  coordinate  system  (which  will  be 
referred  to  as  the  global  system).  The  global  system  chosen  is  two 
translations  and  one  rotation  per  story.  The  origin  of  these  global 
displacement  coordinates  at  each  story  level  is  taken  at  the  center  of 
mass  of  that  story  segment.  This  position  may  vary  from  story  to  story. 
Such  a  formulation  will  degenerate  the  mass  matrix  to  a  diagonal  form,  thus 
simplifying  the  eigen-value  problem  in  the  dynamics. 


The  first  step  is  to  develop  the  transformation  between  the  frame  lateral 
displacements  and  the  global  displacements.  With  reference  to  Figure  11, 
the  transformation  at  any  level,  n,  is  as  follows: 


fr  ^ 
I  xn  i 


r.  =  <cosa  sina  -d  > 
Ln  n 


yn 


v  en 


J 


or: 


r.  =  a  r 
Ln  -n  -n 
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Assembling  the  transformations  for  all  floors,  we  obtain  the  complete 
transformation  between  frame  lateral  displacements  and  global  displace¬ 
ment  as  follows: 


or: 

:u  =  »,  r 

r  is  the  complete  vector  of  global  displacements.  The  frame  lateral 
stiffness  is  transformed  to  the  global  system  and  becomes: 


where  the  subscript  i  denotes  the  ith  frame. 


The  structure  lateral  stiffness  is  assembled  by  the  addition  of  components 
from  all  frames:  i  .e. , 


K  =  I  K. 
i  _1 

The  frame  lateral  load  vector  from  sidesway  effects  must  also  be  trans¬ 
formed  to  the  global  system.  This  transformation  is  shown  by: 
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The  global  load  vector  is  formed  by  the  summation  of  frame  sway  effects 
and  the  addition  of  externally  applied  lateral  loads  F,  i.e.: 

R  »  E  R.  +  F 


The  global  forces  F  are  specified;  however,  they  are  also  given  by: 


*-  -  1  Ai  Pli 


Expanding  the  tri -matrix  product: 
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It  is  worth  noting  that  a  typical  3x3  submatrix  K. ,  within  K  has  the 

1  J  "“1 

form  a.  k.  .  a  ..  Obviously  this  product  may  be  formed  independently  for 
each  term  in  and  added  directly  into  K.  Hence  the  global  equilibrium 
equations  are  formed. 


R  =  K  r 


It  may  be  noted  that  the  global  sviffness  K  is  a  full  matrix,  but  it  is 
of  course  relatively  small  compared  to  the  total  number  of  degrees  of 
freedom  associated  with  all  the  frames  in  the  structure. 
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CHAPTER  III:  STATIC  ANALYSIS 


The  static  analysis  equations: 

R  =  k  r 

are  solved  directly  by  Gaussian  elimination  giving  a  vector  of  global 
lateral  displacements,  r.  Next,  for  each  frame,  the  lateral  displace¬ 
ments,  rLi  are  computed  using: 

rU  .  fli  r 

To  complete  the  solution  for  each  frame,  the  following  system  is  considered. 


Note  that  these  are  the  equations  which  were  reduced,  then  saved  at  each 
level,  n,  of  the  frame.  That  is,  was  triangularized.  At  any  stage, 
n,  rn+.j  and  rL  are  known  and  so  rn  is  computed  by  back  substitution. 

To  start  this  sequence,  we  simply  note  that  for  n  =  N  (the  number  of 
stories  in  the  structure)  r^+^  represents  the  displacements  at  the 
foundation  which  are  zero  since  columns  are  assumed  rigidly  connected  to 
the  foundation.  Thus  the  frame  joint  displacements  are  computed  succes¬ 
sively  story  by  story  and  individual  member  forces  may  be  computed  at 
the  same  time  from  the  force /deformation  transformations  previously 
presented. 
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A.  Vertical  Loads  Analysis 


The  vertical  loads  are  applied  on  each  individual  frame  as  beam  span  loads. 
Four  independent  vertical  loading  conditions  are  possible.  The  self 
weight  of  the  frames  can  be  automatically  calculated  by  the  program  and 
added  to  the  load  vector  of  the  first  load  condition.  Typically  the 
first  load  condition  is  used  for  the  dead  load  analysis  of  the  structure; 
the  second  load  condition  is  used  for  the  live  load  analysis  of  the 
structure.  The  third  and  fourth  load  conditions  may  be  used  for  skip 
live  loading  or  left  unused. 

B.  Lateral  Load  Analysis 

The  lateral  static  loads  are  applied  as  forces  acting  at  a  particular 
point  on  each  floor  level.  Two  independent  lateral  loading  conditions 
are  possible.  The  lateral  loads  may  be  due  to  wind  or  earthquake.  The 
wind  loads  have  to  be  calculated  and  input  by  the  user,  based  upon  the 
wind  pressure  and  the  exposed  tributary  area  of  the  building  at  each 
level  of  the  structure.  The  seismic  static  equivalent  loads  ma^  je 
automatically  calculated  by  the  program,  based  upon  the  requirement  of 
Reference  14.  The  modal  participation  factors  calculated  by  the  program 
are  used  to  determine  the  predominant  directions  of  the  modes  and  the  time 
periods  of  the  predominant  modes  are  used  in  calculating  the  seismic  loads 
in  the  corresponding  directions. 

The  program  has  options  to  calculate  the  dynamic  properties,  such  as  the 
mass  and  mass  moment  of  inertia  of  each  floor  level  based  upon  simplified 
user  input. 
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CHAPTER  IV:  DYNAMIC  ANALYSIS 


The  exact  formulation  of  the  dynamic  response  of  a  structure  involves  an 
infinite  number  of  degrees  of  freedom.  For  most  structures,  however,  the 
response  may  be  adequately  captured  by  a  limited  number  of  discrete 
points  (or  joints)  within  the  system.  In  the  buildings  considered  here, 
the  response  may  be  described  by  the  lateral  motions  of  each  floor  level, 
as  previously  described  for  the  formation  of  the  lateral  structure  stiff¬ 
ness  matrix.  The  center  of  mass  is  used  as  the  master  constraint  location 
at  each  level  in  the  generation  of  the  lateral  stiffness  matrix.  The 
tributary  mass  of  each  story  level  is  lumped  at  the  center  of  mass  of  the 
level  along  with  the  mass  moment  of  inertia  of  the  floor  about  a  vertical 
axis  through  the  center  of  mass  to  compensate  for  the  rotational  aspects 
of  the  lumping  process.  The  resulting  mass  matrix  is  of  diagonal  form. 
With  this  lumped  parameter  idealization,  equilibrium  of  the  structure  is 
described  by  a  set  of  ordinary  second  order  differential  equations. 

A.  Mass  Approximation,  Mass  and  Mass  Moment  of  Inertia 

In  the  diaphragm  shown  in  Figure  12,  there  are  various  lumped  masses 
(m^ ,  m^,  m^  .  .  .  etc.)  and  other  distributed  masses  associated  with 
the  diaphragm  level. 

When  the  diaphragm  is  subjected  to  a  unit  translational  acceleration  in 
the  Y-direction,  inertia  forces  opposing  the  direction  of  the  accelera¬ 
tion  will  be  generated,  i.e.  f^  =  m^  x  1 ,  f^  =  m^  x  1 ,  f^  =  m^  x  1  .  .  . 
The  resultant  of  all  these  forces  and  line  of  action  of  the  resultant 
can  be  determined.  The  magnitude  of  the  resultant  is  found  to  be  =  f-j  + 
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*2  +  *3  ‘  ‘  *  =  +  m2  +  m3  1  '  ‘  =  tota^  mass  associated  with  the 

diaphragm.  The  resultant  is  parallel  to  the  V-direction  and  passes 
through  a  point  at  a  distance  Xm  from  0. 

Similarly,  a  unit  translation  in  the  X-direction  will  give  a  resultant 
of  the  same  magnitude  but  parallel  to  the  X-direction  and  passing 
through  a  point  a  distance  Ym  from  0. 

The  coordinates  X  ,  Y  define  the  location  of  a  point  known  as  the 
m  m  r 

center  of  mass. 

Redefining  the  term.  Mass:  "The  mass  of  a  diaphragm  nay  be  defined  as 
the  force  generated  when  the  center  of  mass  of  the  diaphragm  under¬ 
goes  a  unit  translational  acceleration.  This  force  acts  at  the  center 
of  mass,  resulting  in  no  associated  moment." 

n 

Mass  =  z  m. 
i=1  1 

Similarly,  defining  the  term.  Mass  Moment  of  Inertia  (or  Rotational 
Mass):  "The  mass  moment  of  inertia  of  a  diaphragm  may  be  defined  as 
the  moment  generated  when  the  center  of  mass  of  the  diaphragm  under¬ 
goes  a  unit  rotational  acceleration  about  a  vertical  axis.  No  resultant 
translational  force  is  associated  with  the  couple." 

The  radial  distances  from  the  center  of  mass  of  the  lumped  masses 

m^ ,  m,,,  m3  .  .  .  are  d^ ,  d£,  d^  .  .  .  respecti vely,  as  shown  in  Figure  12. 

Due  to  a  unit  rotational  acceleration  of  the  center  of  mass  about  a 
vertical  axis,  ,  m^,  .  .  .will  have  translational  accelerations 

of  d^xl,  d^xl ,  d^xl  .  .  .  Thereby  giving  corresponding  inertia  forces  of 
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m-|d^,  m^d^,  m^d^  ....  The  moments  of  these  forces  about  a  vertical 

2  2  2 

axis  through  the  center  of  mass  are  m^d^ ,  11^2,  m^  .... 

MMI  =  m^d^  +  n^d^  +  m^d^  .... 

n  2 
=  r  m.d. 
i=l  1  1 

=  Polar  Moment  of  Inertia  of  all  Masses,  about  a 
vertical  axis  through  the  center  of  mass 


B.  Dynamic  Equilibrium  Equations 


The  equilibrium  equations  for  a  structure,  including  dynamic  effects, 
may  be  written  in  the  following  form: 


Hra  +  C  r  +  K  r  =  P  (t) 


(a) 


where:  M 
C 
K 

P  (t) 
r 


mass  matrix 
damping  matrix 
stiffness  matrix 

applied  load  vector,  which  may  be  time  dependent 

displacement  vector  of  deformation  relative  to 
support  motion 

absolute  acceleration  vector 
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r  and  ra  are  related  in  the  following  fashion: 


r  =  v„  +  r 
-a  -g 


where  is  the  vector  of  pseudo-static  displacements  due  to  support 
movement.  Also: 


r  =  v  +  r 
-a  -g 

These  vectors  have  the  following  form  for  a  typical  floor,  of  a  building 
shown  in  Figure  13. 
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Or,  for  all  floors: 


r  =  B  V  +  r 
-a  -  g 


er: 


X 
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Figure  13.  Ground  and  structural  displacements 
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where : 


/  \ 

*1 

-2 

B  =  \  b3  >  ;  etc. 

^  -N  > 

In  the  case  of  seismic  analysis,  there  are  no  externally  applied  loads; 
i.e.,  P  (t)  =  0.  Then  equation  (a)  may  be  written  as; 

M  (r  +  B  v  )  +  C  r  +  Kr  =  0 

-  -  -  9  . 

or: 

Mr  +  Cr  +  Kr=-MBv'g  . (b) 

This  coupled  set  of  equations  may  be  solved  simultaneously  with  an 
appropriate  numerical  technique.  Another  approach,  which  will  be  used 
here,  is  to  find  a  transformation  which  uncouples  the  equations  so  that 
they  may  be  solved  independently.  This  transformation,  of  course  is  via 
the  eigen-vectors  or  mode  shapes  of  the  system. 

C.  Mode  Shapes  and  Frequencies 

The  vibration  mode  shapes  represent  the  solution  of  the  undamped  free 
vibration  problem  given  by 

Mf  +  Kr  =  0 
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The  eigen-value  problem  to  be  solved  is  written  as: 


K  $ 


whert-  $  =  mode  shapes 
ui  =  frequencies 


The  mode  shapes  are  normalized  such  that: 

^  H  ^  =  I 


then  also: 


^  K  | 


2 

UI 


Also,  it  is  assumed  that  the  damping  matrix  C  is  of  a  form  that  is 
uncoupled  by  the  mode  shapes;  specifically  it  is  assumed  that: 

<(>  C  4  =  [  2  X  u>  ] 

so  that  represents  the  damping  of  the  m  th  mode. 


The  actual  displacements,  r,  are  now  expressed  as  a  linear  combination 
of  the  mode  shapes. 


r  = 


[  $1  ^2  -3  '  '  -N  1 


z^tr 

Zo(t) 


ZN^- 


•  •  (c) 


i .e. :  r  =  $  Z 

also  r  =  £  1 
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and : 


r  =  }  Z 

where  Z  (t)  represents  the  response  of  the  m  th  mode. 


D.  Time  History  Analysis 

Using  equations  (  b  )  ,  equation  (  c  )  may  be  rewritten  as: 

M  $  Z  +  C  $  2  +  KfZ  =  -  M  B  v 


Premultipl  ication  by  <)J  yields  the  uncoupled  set  of  second  order 
equations : 

M*  1  +  C*  l  +  K*  Z  =  P*  v„ 


where : 


m*  =  m  =  i 

c*  *  /  c  f  ■  [  2>m  ^  ] 

K*  =  *T  K  i  =  [  J  ] 


P*  v 


M  B  v 
-  -  9 


to  find  the  form  of  P*,  consider: 


M  B  = 


*T 

f  \ 

ml 

sinn 

ml 

COSB 

J1 

0 

1^2 

J 

sing 

nip 

\ 

cosg 

0 

0 

V  s 
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where:  m,  =  mass  of  story  1 


=  rotational  mass  moment  of  inertia  of  story  1 


m-j  sins 


m-j  cose 


MB  = 


m2  sine 


m2  cose  i 


:  j 


So,  a  typical  term  of  P*  has  the  form: 


P*  =  4>  ‘  M  B 

m  ^-m 


=  **lx  *ly  *10  *2x  *2y  *28  ' 


>  1  sm3 

■  m^  coss 


m2  sing 


m2  cos  S 


P*  =  l  m  {sins 
"  n=l  " 


>„*  *  cose  1 


1  % 


Now  a  typical  equation  governing  the  response  in  the  m  th  mode  has  the 
form: 


2  X  w  1  + 

m  m  m 


a j2  Z 
m  m 


P*  v 
m  g 


(d) 


For  any  earthquake,  the  ground  acceleration,  vg  is  specified  as  a  set 
of  discrete  values  and  linear  interpolation  is  used  for  intermediate 
values.  On  any  linear  portion  then: 


vg  =  A  +  Bt 


where  A  and  B  are  computed  from  the  end  values  as  shown  in  Figure  14. 


On  any  linear  segment  t-j ,  t^  then: 


Z  +  2\  u  l  +  uf  1  =  P*  (A  +  Bt) 

m  m  m  m  mm  m 


The  solution  to  this  equation  is  summarized  in  Figure  14. 


At  rest  initial  conditions  are  used  for  the  first  linear  portion.  The 
values  of  displacement  and  velocity  at  the  end  of  any  linear  portion 
form  the  initial  conditions  for  the  following  linear  segment  and  so  on. 
Repetition  gives  the  complete  solution  over  the  required  time  span. 

With  solutions  for  each  mode,  equation  (c)  is  used  to  give  a  set  of 
structure  displacements  r  at  each  output  time  step. 

The  backsubstitution  procedure  used  for  the  time  history  analysis  is 
exactly  the  same  as  that  described  for  the  static  analysis  in  Chapter  III. 
Backsubstitution  for  each  time  step  is  equivalent  to  one  static  load 
backsubstitution.  The  frame  displacements  and  member  forces  are 
determined  at  each  time  step  and  the  maxima  of  these  parameters  over  the 
time  span  are  output  as  dynamic  load  condition  3. 
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E.  Response  Spectrum  Analysis 

Unless  actual  histories  of  displacements  and  forces  are  required  for  a 
specific  earthquake  a  more  realistic  and  economical  approach  for  dynamic 
analysis  is  via  the  response  spectrum  method.  For  a  particular  ground 
motion  history  v  (t),  the  spectrum  curve  is  defined  as  follows: 

The  response  of  a  single  mass  system  with  damping  X,  _.nd  circular 
frequency  u>,  subjected  to  a  ground  motion  history  v  (t)  is  governed  by 
the  equation: 

u(t)  +  2Aw  u(t)  +  u)^u(t)  =  Vg ( t) 

Let  u  be  the  maximum  absolute  value  that  u(t)  attains.  A  plot 
of  this  maximum  displacement  versus  the  frequency  w  for  each  A  is 
by  definition  the  displacement  response  spectrum  (Sd)  for  the  earth¬ 
quake  Vg(t).  A  plot  of  is  the  pseudo-veloci ty  spectrum  (PSv) 

and  a  plot  of  u,  u  is  the  pseudo-acceleration  spectrum  (PSa). 

ITlaX 

These  pseudo-velocity  and  acceleration  spectra  are  of  the  same  physical 
interest  but  are  not  an  essential  pa;  t  of  a  response  spectrum  analysis. 

Recalling  equation  (d),  if  the  dynamic  loading  on  the  s 1 1  jre  is 

specified  in  terms  of  the  pseudo-acceleration  spectrum,  then  the  maximum 
response  for  the  m  th  mode  is  given  by: 


Z 

m 


max 


P* 

m 


PSa  (uj  ,  x  ) 
'  m  nr 

U> 
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Therefore,  the  maximum  contribution  of  mode  m  to  the  total  three 
dimensional  response  of  the  structure  is: 

r  =  Z  it 
-m  m  -m 
max 

For  all  modes  Sd  is,  by  definition,  positive.  The  maximum  modal 
displacement  rm  is  proportional  to  the  mode  shape  4>;  and  the  sign 
of  the  proportionality  constant  is  given  by  the  sign  of  the  modal 
participation  factor,  P*.  Therefore,  each  maximum  modal  displa  ient 
has  a  unique  sign.  Also,  the  maximum  internal  modal  forces,  which  are 
consistently  evaluated  from  the  maximum  modal  displacements,  have 
unique  signs. 

A  complete  analysis  is  performed  down  to  the  member  force  level  with 
the  maximum  modal  displacements  of  the  structure  for  each  individual 
mode  using  the  backsubsti tution  procedure  described  in  Chapter  11. 

The  maxima  in  each  mode  will  generally  occur  at  different  times.  The 
combination  of  the  modal  components  of  the  displacements  and  member 
forces  to  give  resultant  values  for  design  purposes  is  performed  at 
the  design  parameter  level  by  the  following  methods. 

1.  The  Square-Root-of-the-Sum-of-the-Squares  (SRSS)^1 2 3^  method 

2.  The  Absolute  Sum  (ABS)^3^  method 

3.  The  Complete  Quadratic  Combination  (CQC)^^  method 

The  SRSS  method  and  the  ABS  method  entirely  neglect  the  signs  of  the 
modal  contributions.  The  SRSS  method  in  general  gives  good  approxi- 
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mations  of  the  dynamic  response  in  structures  with  well  separated 
frequencies.  The  ABS  method  is  basically  for  interest  to  give  an 
upper  bound  on  the  maximum  values. 

In  structures  with  closely  spaced  modes  or  multiple  frequencies,  the 
fact  that  the  SRSS  method  neglects  the  signs  of  the  modal  components 
may  cause  the  design  parameters  to  be  dramatically  overestimated  in 
some  elements  while  being  significantly  underestimated  in  other  elements. 
The  CQC  method  overcomes  this  difficulty  and  it  is  recommended  as  the 
best  of  the  three  methods  for  obtaining  the  most  realistic  results. 

F.  Dynamic  Options 

The  dynamic  options  currently  available  in  CTABS80  are: 

1.  Calculation  of  mode  shapes  and  periods  (frequencies) 

2.  Response  spectrum  analysis  for  any  acceleration  spectrum 
supplied  by  the  user  using  the: 

a.  SRSS  modal  combination  as  Dynamic  load  condition  1 

b.  Sum  of  absolute  value  modal  combinations  as  Dynamic  load 
condition  2 

c.  Complete  Quadratic  combinations  as  Dynamic  load  condition  3 

3.  Time  history  analysis  maxima  for  any  ground  motion  supplied  b} 
the  user  as  Dynamic  load  condition  3 

Either  dynamic  analysis  condition  may  be  combined  with  any  static  load 
condition. 
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1.  SRSS  Combination 


(Ixl)  -^/(Ixn)  (nxn)  (nxl) 


Where  is  an  identity  matrix 


2.  ABS  Combination 


F  =  fT  sign  f 
(1x1)  (Ixn)  (nxl) 


Where  sign  f  is  a  unit  matrix  contain¬ 
ing  the  signs  of  the  corresponding 
elements  of  matrix  f 


3.  CQC  Combination 


F 

(lxl) 


C 

(nxn) 


f 

(nxl) 


NOTES/ 

f  =  vector  of  modal  components 
F  =  combined  resultant 
n  =  number  of  modes 


Where  C  is  the  matrix  of  modal  cross¬ 
correlation  coefficients  given  by: 

C  =  8X2(HT)r3/7 

ij  (l-r2)2+4x^r(l+r)2 

where  r  =  the  ratio  of  the 

circular  frequencies  of  the  coupling 

modes  and  x  is  the  damping  associated 

with  the  response  spectrum  curve 

being  used. 


Figure  15.  Summary  of  modal  combination  techniques  used  in  CTABS80 
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CHAPTER  V:  GENERAL  OBSERVATIONS 


A.  Program  Appl ication 

The  effective  application  of  a  computer  program  for  the  analysis  of 
practical  situations  involves  a  considerable  amount  of  experience.  The 
most  difficult  phase  of  the  analysis  is  in  assembling  an  appropriate 
model  which  captures  the  major  characteristics  of  the  structural 
behavior  of  the  building.  No  computer  program  can  replace  the  engineering 
judgement  of  an  experienced  engineer.  It  is  well  said  that  an  incapable 
engineer  cannot  do  with  a  ton  of  computer  output  what  a  good  engineer  can 
do  on  the  back  of  an  envelope.  Correct  output  interpretation  is  just  as 
important  as  the  preparation  of  a  good  structural  model.  Verification  of 
unexpected  results  needs  a  good  understanding  of  the  basic  assumptions  and 
the  mechanics  of  the  program.  Static  equilibrium  checks  are  necessary  not 
only  to  check  the  computer  output  but  to  understand  the  basic  structural 
behavior  of  the  building. 

B.  Static  Seismic  Analysis  of  Buildings 

At  the  present  time,  the  seismic  design  of  most  buildings  in  California 
and  other  earthquake  regions  of  the  United  States  is  based  upon  the 
Uniform  Building  Code.  The  UBC  method  allows  the  seismic  loads  to  be 
approximated  by  an  equivalent  set  of  lateral  static  loads.  The  magnitude 
of  the  loads  is  based  upon  the  seismic  zone,  the  structural  system,  and 
the  fundamental  period  of  the  structure.  Corrections  to  compensate  for 
local  soil  conditions  and  the  physical  importance  of  the  structure  are 
also  defined. 
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An  approximate  formula,  specified  in  the  UBC,  may  be  used  to  estimate  the 
fundamental  period.  The  period  associated  with  the  predominant  structural 
mode  obtained  via  the  TABS  program  is  more  accurate  and  appropriate.  The 
suggested  UBC  distribution  of  the  lateral  loads  over  the  height  of  the 
building  is  triangular  with  some  correction  to  allow  for  higher  mode 
effects.  Behavior  of  structures  that  have  dynamically  decoupled  regions 
due  to  stiffness  and/or  mass  discontinuities,  causing  significantly  non- 
tri angular  inertia  load  patterns  are  not  adequately  covered  by  the  code. 

By  examining  the  structural  modes  produced  by  a  TABS  analysis  such 
structural  complexities  can  be  isolated. 

The  determination  of  the  minimum  horizontal  torsional  design  moments,  as 
specified  by  the  UBC  for  the  design  of  structures  having  rigid  diaphragms, 
requires  the  location  of  a  center  of  rigidity  of  the  structure  at  each 
level.  The  definition  of  torsional  moments  on  such  a  basis  for  multi¬ 
story  structures  is  technically  vague.  It  is  only  meaningful  in  single 
story  structures  where  there  are  no  stories  above  or  below  to  affect  the 
rotation  of  the  level  under  consideration. 

The  UBC  lateral  loads  are  only  a  small  fraction  of  the  loads  developed 
during  a  significant  earthquake,  and  must  therefore  be  considered  as 
minimum  requirements.  As  a  result  of  the  above  mentioned  inadequacies 

the  need  for  a  more  comprehensive  code  earthquake  analysis  methodology  is 

(8) 

apparent  to  most  structural  engineers'  . 

C .  Computer  Methods  Versus  Hand  Methods 

High  speed  digital  computers  and  the  development  of  computer  programs  such 
as  TABS  have  given  engineers  the  capability  to  consider  aspects  of 
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structural  behavior  that  conventional  hand  analysis  techniques  have 
traditionally  neglected.  Hand  analysis  techniques  used  by  practicing 
engineers  for  the  lateral  analysis  of  multistory  structures  have  been 
shown  to  violate  joint  statics  and  compatibility. 

The  following  examples  demonstrate  the  degree  of  error  that  could  be 
present  in  a  conventional  hand  analysis,  by  comparison  with  a  TABS 
analysis;  i.e.,  one  that  completely  satisfies  statics,  compatibility, 
and  boundary.  Examples  presented  are  of  simple  symmetrical  multistory 
buildings  with  symmetrical  loading.  The  stories,  therefore,  translate 
under  lateral  load  without  rotating,  thereby  keeping  the  problems  clear 
and  demonstrative.  The  proportions  of  the  structures  and  magnitudes  of 
the  forces  have  been  chosen  to  generate  problems  of  a  nature  that  a 
conventional  structural  engineer  is  commonly  faced  with  in  practice. 

(i).  Example  1 

This  is  a  classic  example  of  shear  wall-frame  interaction.  A  10  story 
shear  wall  is  connected  in  parallel  with  a  ten  story  frame  at  each  story 
level  through  a  rigid  link.  The  axial  deformations  in  the  beams  are 
neglected,  thus  simulating  a  rigid  diaphragm.  Therefore,  the  lateral 
displacements  of  the  respective  stories  of  the  frame  and  shear  wall  are 
equal.  See  Figure  16. 

Consider  the  top  story.  Based  on  a  conventional  hand  analysis,  one 
would,  in  general,  tend  to  ignore  the  stiffness  of  the  frame  and  conclude 
that  the  shear  wall  takes  close  to  100%  of  the  applied  100  K,  and  that  the 
frame  being  relatively  flexible  gets  a  negligible  amount  of  the  shear. 
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A  "correct"  analysis,  however,  reveals  that  the  frame  has  a  total  shear 
of  186  K  at  that  level.  Thus,  the  frame,  besides  carrying  the  total  100  K 
applied  load,  is  laterally  supporting  the  shear  wall  which  puts  an  addi¬ 
tional  86  K  on  the  frame.  The  shear  wall  is  in  effect  acting  as  a  propped 
cantilever  supported  by  the  frame  at  the  upper  story  levels. 

The  phenomenon  may  be  explained  as  follows:  If  the  frame  and  shear  wall 
are  loaded  independently  with  the  load,  and  the  lateral  displacements  of 
the  respective  floors  compared,  it  will  be  observed  that  the  shear  wall 
has  larger  lateral  displacements  in  the  upper  levels,  whereas  the  frame 
has  larger  displacements  in  the  lower  stories  and  vice  versa.  Compati¬ 
bility  of  joint  rotations  will  have  a  significant  effect  on  these  dis¬ 
placement  patterns.  When  loaded  together,  the  constraint  of  equal  story 
displacements  is  enforced,  thus  resulting  in  this  unique  shear  distribu¬ 
tion.  Notice  that  in  the  lower  stories  the  shear  gradually  shifts  to 
the  shear  wall. 

This  example  demonstrates  the  importance  of  the  interaction  of  all  the 
elements  on  one  another  and  that  the  hand  analysis  method  of  analyzing 
an  n-story  structure  as  n  1 -story  structures  stacked  one  over  the  other 
with  no  interaction  of  one  on  the  other  can  lead  to  highly  unreliable 
results. 

(ii ) .  Example  2 

This  problem  consists  of  two  2-story  walls  and  one  1 -story  wall.  See 
Figure  17a.  Again,  as  in  Example  1,  the  walls  are  connected  by  rigid 
links  at  the  floor  levels  to  simulate  a  rigid  diaphragm  and  to  enforce 
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equal  lateral  story  displacements  in  all  the  walls.  The  walls  are  loaded 
with  a  total  story  load  of  50  K  at  each  level.  As  expected  in  a  convention¬ 
al  hand  analysis,  the  2-story  walls  equally  share  the  50  K  lateral  load  at 
the  upper  level . 

Now  let  us  consider  the  shear  distribution  in  the  bottom  story.  The  base 
shear  is  100  K.  In  the  light  of  the  fact  that  there  are  three  resisting 
elements  at  this  level,  all  40'  long  and  1'  thick,  a  conventional  hand 
analysis  would  conclude  that  the  base  shear  will  be  carried  equally  by  all 
three  elements;  that  is  33  K  each. 

A  "correct"  analysis,  however,  indicates  that  the  1-story  wall  takes  over 
50%  more  shear  than  each  2-story  wall.  In  Figure  17b  are  presented  the 
free  body  diagrams  of  the  lower  levels  of  the  1-story  wall  and 
the  2-story  walls.  Consider  the  lateral  story  displacements  of  the  1st 
level  in  each  wall.  In  the  1-story  wall  the  lateral  displacement,  aP ^  ,  is 
due  to  ,  the  shear  force  in  the  wall.  In  the  2-story  wall  the  lateral 
displacement  is  due  to  two  factors.  Firstly,  AP^,  that  is  due  to  P^.  the 
shear  force  in  the  wall  and,  secondly,  AM2 ,  due  to  M^,  the  moment  at  the 
top  of  the  wall  due  to  the  fact  that  the  wall  is  2  stories  high.  Now  for 
the  lateral  displacements  to  be  equal: 

AP-|  =  AP2  +  AM2 

Therefore  aP]  >  aP2 

so  that  P.j  >  P2 

The  discrepancy  between  the  conventional  hand  analysis  method  and  the 
"correct"  method  here,  again,  is  due  to  the  fact  that  the  effect  of  the 
upper  story  on  the  lower  story  is  accounted  for  incompletely. 
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( i i i ) .  Example  3 


This  structure  consists  of  4  walls.  Two  4-story  walls  and  two  2-story 
walls.  See  Figure  18a.  Again,  the  walls  are  connected  by  rigid  links  at 
the  floor  levels  to  simulate  a  rigid  diaphragm. 

In  the  "correct"  analysis  the  shear  forces  in  the  10  foot  walls  in  the 
3rd  and  4th  levels  are  as  would  be  expected  in  a  conventional  hand 
analysis.  Note  the  shear  distribution  at  the  2nd  level.  At  this  level 
there  is  a  considerable  increase  in  the  story  stiffness  due  to  the  two 
40-foot  walls.  This  restricts  the  lateral  diaphragm  movement  to  the 
extent  that  the  10‘  walls  are  in  effect  laterally  supported  by  the  dia¬ 
phragm  at  this  level,  and,  therefore,  behave  like  over-hanging  cantilevers 
as  shown  in  Figure  18b.  This  explains  why  these  walls  have  a  negative 
shear  of  154  K  each  at  this  level. 

The  conventional  hand  analysis  method  for  such  problems  completely  dis¬ 
regards  the  possibility  of  negative  shear  forces  occurring  in  the  walls, 
thereby  always  assuming  that  the  walls  support  the  diaphragm  laterally, 
and  not  recognizing  that  at  times  the  diaphragm  may  actually  be  the 
support  for  certain  walls  at  certain  levels. 

(iv) .  Example  4 

This  example  demonstrates  the  effect  of  axial  deformations  on  the 
distribution  of  shear  to  a  series  of  walls.  Consider  the  structure 
shown  in  Figures  19a and  19b.  The  structure  consists  of  two  solid  walls, 
and  one  wall  terminating  on  two  columns.  In  case  A,  the  columns  are 
12"  square,  see  Figure  19a.  Again,  the  walls  are  connected  by  rigid 


links  at  the  story  levels. 


Consider  the  shear  distribution  in  the  upper  level.  In  case  A,  a 
"correct"  analysis  indicates  that  the  50  K  story  shear  is  distributed 
approximately  equally  among  the  3  walls.  A  conventional  hand  analysis 
would  lead  us  to  the  same  conclusion.  In  case  B,  however,  the  shear 
taken  by  the  wall  terminating  on  columns  is  50%  of  that  taken  by  the  solid 
walls.  The  reason  being  that  the  smaller  axial  area  of  the  columns  gives 
larger  axial  deformations  which,  in  turn,  reduce  the  lateral  rigidity  of 
the  wall  in  the  story  above. 

A  conventional  hand  analysis  neglects  the  effects  of  axial  deformations , 
and,  therefore,  would  give  an  equal  shear  distribution  in  all  three  walls 
regardless  of  the  size  of  the  columns. 

(v) .  Example  5 

This  example  clearly  demonstrates  the  analytical  discrepancy  in  analyzing 
an  n-story  structure  as  n  1 -story  structures.  The  structure  is  a  shear 
wall  with  the  same  story  height,  wall  dimensions  and  openings  at  every 
level .  See  Figure  20. 

A  conventional  hand  analysis  (analysis  of  the  structure  as  four  1-story 
structures)  would  indicate  that  the  location  of  the  point  of  contra- 
flexure  of  the  piers  in  all  the  stories  would  be  at  a  constant  distance 
vertically  from  the  corresponding  diaphragm.  Also,  the  ratio  of  the 
shear  force  in  pier  A  to  that  in  pier  B  in  all  stories  will  be  constant. 

A  "correct"  solution  of  the  structure,  however,  shows  that  there  is  no 
point  of  contra-flexure  in  pier  A  in  the  first  three  stories.  Also,  the 
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percentage  of  the  story  shear  carried  by  pier  B  is  not  constant  in  all 
stories  but  decreases  as  we  move  into  the  lower  stories. 

Joint  rotation  incompatibility  from  story  to  story  is  the  main  cause  for 
this  discrepancy  between  the  hand  analysis  and  the  "correct"  analysis. 

( vi ) .  Example  6 

This  is  another  example  demonstrating  the  effect  of  axial  deformations 
on  the  shear  distribution.  The  structure  consists  of  slender  piers 
framing  into  relatively  stiff  spandrels  and  is  8  stories  tall.  See 
Figure  21.  For  all  practical  purposes,  the  piers  may  be  considered  fixed 
in  rotation  at  both  ends.  Since  all  the  piers  are  of  the  same  size,  a 

conventional  hand  analysis  would  indicate  that  the  story  shear  is 

distributed  equally  among  the  5  piers. 

A  "correct"  analysis,  as  we  can  see,  indicates  that  the  piers  closer  to 
the  center  take  a  higher  percentage  of  the  story  shear.  In  the  top 
story  for  instance,  the  center  pier  takes  over  65%  higher  shear  than  the 

end  pier.  This  discrepancy  is  due  to  the  fact  that  the  "correct" 

analysis  considers  axial  deformations  in  the  piers,  whereas  the  hand 
analysis  does  not. 

This  behavior  may  be  explained  by  the  following  analogy.  Consider  the 
lateral  displacements  in  a  vertical  cantilever  with  a  rectangular  cross 
section  and  lateral  loading.  If  the  shear  deformations  are  negligible 
compared  to  the  bending  deformations  in  the  cantilever,  the  distribution 
of  the  shear  stress  across  the  section  is  parabolic  with  the  maximum  at 
the  center.  However,  if  the  deformation  pattern  is  one  of  pure  shear. 
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such  that  the  bending  deformations  are  negligible,  the  shear  stress  dis¬ 
tribution  is  constant  over  the  section  and  not  parabolic. 


The  axial  deformations  of  the  piers  in  this  example  correspond  to  the 
bending  deformations  in  the  analogy  above.  If  the  piers  were  axially 
stiff  to  the  extent  that  there  were  no  axial  deformations,  the  shear 
distributions  in  the  piers  would  be  equal.  However,  as  the  piers  deform 
axially,  the  piers  closer  to  the  middle  have  heavier  shear.  This  also 
explains  why  the  shear  distribution  is  more  uniform  as  we  move  into  the 
lower  stories. 

D .  Dynamic  Seismic  Analysis  of  Buildings 

The  deficiencies  of  the  present  seismic  design  procedures  are  clearly 
summarized  in  Reference  8.  It  is  apparent  that  the  present  code  is  a 
very  approximate  method  based  on  the  first  mode  only.  The  foundation 
factors  discussed  later  are  not  considered.  Another  factor  which  is 
important  in  an  elastic  analysis  is  the  damping.  Spectra  for  damping 
of  2  and  10*  are  shown  in  Figure  22.  It  is  clear  that  the  Uniform 
Building  Code  seismic  loads  are  very  small  compared  to  the  force.,  pro¬ 
duced  in  recorded  earthquakes.  It  has  been  estimated  that  earthquakes 
of  the  Parkfield  magnitude  can  be  expected  about  once  per  year  at  some 
point  in  California,  and  earthquakes  of  the  El  Centro  magnitude  may  be 
expected  every  five  or  six  years. 

The  selection  of  a  design  spectrum  for  the  response  spectrum  analysis  of 
a  particular  building  will  depend  on  the  geographical  area,  the  local  soil 
condition,  the  type  of  construction  material  and  the  intended  use  of  the 
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Figure  22.  Response  spectra 


building.  Many  Soils  Engineering  firms  now  specialize  in  the  dynamics  of 
soil  systems  to  evaluate  specific  sites  and  recommend  shapes  and  intensi¬ 
ties  for  dynamic  response  spectra.  Most  lending  agencies  are  requiring 
dynamic  analysis  of  structures  as  part  of  their  financing  terms  for 
buildings  in  major  metropolitan  areas  of  California. 

For  certain  types  of  earthquakes  it  has  been  observed  that  the  vertical 
accelerations  are  comparable  in  magnitude  to  the  lateral  accelerations . 
However,  all  building  have  been  designed  elastically  for  a  minimum  of 
1  g  in  thevertical  direction;  therefore,  these  additional  vertical  forces 
very  often  do  not  cause  direct  damage  to  the  structure.  Of  course,  they 
should  be  considered  in  the  design  of  members  in  addition  to  the  lateral 
earthquake  loads.  For  most  structures  the  stiffness  in  the  vertical 
direction  is  very  large;  hence,  the  vertical  periods  will  be  very  small. 
Therefore,  a  dynamic  analysis  in  the  vertical  direction  may  not  be 
required.  A  direct  increase  in  dead  load  stresses  may  be  a  good  method 
to  approximate  the  effects  of  vertical  earthquake  loads. 

E.  Foundation  -  Building  Interaction 

Foundation  modeling  has  always  been  an  area  of  particular  concern.  The 
vertical  and  rotational  stiffnesses  under  each  column  can  be  easily  input 
by  providing  an  extra  "dummy"  story.  However,  the  assigning  of  accurate 
stiffness  values  for  these  soil  springs  can  be  difficult. 

In  recent  years  considerable  research  has  been  conducted  in  the  area  of 
foundation  -  building  interaction.  However,  very  little  of  this  work  has 
been  of  direct  value  to  the  profession  involved  in  the  earthquake  analysis 


63 


of  buildings.  Several  of  the  suggested  approaches  have  been  difficult  to 
apply  in  case  of  complex  buildings,  or  they  have  had  serious  theoretical 
restrictions. 

Before  foundation  interaction  effects  are  included  in  the  analysis  it  is 
necessary  to  define  the  exact  location  of  the  earthquake  input.  If  the 
design  criteria  states  that  the  input  is  at  the  base  of  the  building  then 
it  is  impossible  to  say  that  the  building  will  modify  the  input,  and  it 
is  impossible  to  include  interaction  effects. 

A  large  amount  of  research  in  this  area  has  been  associated  with  machines 
vibrating  on  an  infinite  foundation  where  the  term  radiation  damping  has 
been  used.  This  work  has  little  value  in  earthquake  engineering  since  the 
energy  source  is  not  at  the  base  of  the  building.  It  is  easy  to  show  that 
the  energy  stored  in  the  building  is  very  small  compared  to  the  energy 
stored  in  the  immediate  foundation  area  in  the  case  of  earthquake  input. 
Also,  the  machine  vibration  problem  is  a  steady  state  phenomenon,  whereas 
earthquakes  produce  a  transient  loading. 

The  continuous  foundation  contains  an  infinite  number  of  degrees  of 
freedom.  Therefore,  any  approach  which  suggest  representing  the  lateral 
behavior  of  the  foundation  with  a  simple  spring,  dashpot  and  mass  system 
is  a  very  gross  approximation.  In  fact,  this  technique  can  produce  a 
filtering  effect  on  the  earthquake  input  and  cause  serious  errors.  For 
lateral  earthquake  input,  this  type  of  approximation  is  only  acceptable 
in  the  representation  of  the  rotational  stiffness  at  the  base  of  columns 
and  shear  walls. 


The  most  significant  factor  to  consider  is  the  modification  of  the  basic 


CHAPTER  VI:  INTERNAL  ORGANIZATION 

An  outline  of  the  subroutine  structure  of  CTABS80  is  presented  in 

Figure  24.  There  are  seven  major  calls  from  CTABS80  associated  with 

the  seven  major  blocks  of  the  program. 

1.  The  first  operation  is  to  read  the  basic  control  information.  The 
data  associated  with  the  complete  building  (story  data  and  structural 
lateral  loads)  is  then  rolled  in  via  subroutine  TABI. 

2.  The  next  operation  involves  reading  in  the  frame  data  of  every  dif¬ 
ferent  frame  in  the  structure.  The  frame  elevations  are  plotted,  if 
requested.  In  non-data  check  modes  the  frame  stiffnesses  are  formu¬ 
lated  and  reduced  and  the  frame  lateral  stiffness  matrices  and  back- 
substitution  equations  are  written  sequentially  on  disc.  This 
operation  is  implemented  by  the  call  to  subroutine  TABF. 

3.  The  call  to  subroutine  TABL  reads  the  frame  location  data  and  formu¬ 
lates  the  complete  lateral  structural  stiffness  matrix  of  the  whole 
building. 

4.  Subroutine  SFRAME  causes  a  plan  view  of  the  building  to  be  plotted, 
showing  the  frame  locations  and  the  directions  of  their  local  axes. 

5.  The  call  to  subroutine  TABE  gives  the  modeshapes  and  freouencies 
of  the  structure  (TABM)  and  triggers  the  automatic  UBC  lateral 
seismic  load  calculation  (TUBC) .  Also  the  dynamic  analysis  control 
information  is  read  in  by  this  call.  Structural  lateral  displace¬ 
ments  due  to  the  static  loads  (TABQ)  and  response  spectrum  dynamic 
loads  are  obtained  at  this  stage. 
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6.  Subroutine  TABDY  reads  in  the  time  history  earthquake  ground  motion 
data  and  causes  the  structural  lateral  displacements  to  be  calculated 
for  each  time  step. 

7.  Finally,  subroutine  TABP  is  called.  This  subroutine  calls  TABC  to 
read  the  load  case  definition  data  for  each  frame.  Then  TABU  is 
called  to  print  the  frame  lateral  displacements  and  TABO  is  called 
to  calculate  the  frame  joint  displacements  for  each  static  load 
condition  and  each  spectral  mode  or  response  time  increment 

from  the  back  substitution  equations  previously  saved. 

As  the  displacements  are  calculated  the  member  forces  are  also 
evaluated  and  printed  according  to  the  load  case  definition 
data  (TABW). 
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CHAPTER  VII:  CONCLUSIONS 


A  general  computer  program  for  the  elastic  three-dimensional  static  and 
dynamic  analysis  of  frame  and  shear  wall  buildings  has  been  presented. 

For  buildings  which  can  be  approximated  by  independent  frames  and  shear 
walls  the  program  is  very  economical  and  easy  to  use  as  compared  to  a 
general  purpose  three-dimensional  structural  analysis  program. 

Many  new  options  have  been  implemented  in  this  release  to  make  the  pro¬ 
gram  a  more  practical  and  useful  engineering  tool. 

The  program  is  based  on  linear  theory.  Non-linear  behavior  such  as  P-A 
effects  and  material  plasticity  are  not  captured  by  the  program. 

If  non-linear  effects  are  to  be  considered  a  step-by-s^ep  response 
analysis  is  required;  however,  this  involves  a  significant  increase  in 
computational  effort  and  will  be  justified  for  only  a  limited  number  of 
buildings.  In  addition,  the  non-linear  material  properties  both  for  most 
structural  and  non-structural  members  have  not  been  established  accurately 
from  experimental  work. 
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